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Two objects can be distinguished if they have different measurable properties. Thus, distinguishabUity 
depends on the Physics of the objects. In considering graphs, we revisit the Ising model as a framework to 
define physically meaningful spectral invariants. In this context, we introduce a family of refinements of the 
classical spectrum and consider the quantum partition function. We demonstrate that the energy spectrum 
of the quantum Ising Hamiltonian is a stronger invariant than the classical one without refinements. For the 
purpose of implementing the related physical systems, we perform experiments on a programmable 
annealer with superconducting flux technology. Departing from the paradigm of adiabatic computation, we 
take advantage of a noisy evolution of the device to generate statistics of low energy states. The graphs 
considered in the experiments have the same classical partition functions, but different quantum spectra. 
The data obtained from the annealer distinguish non-isomorphic graphs via information contained in the 
classical refinements of the functions but not via the differences in the quantum spectra. 



Kac's' question "Can one hear the shape of a drum?" is part of the scientific pop culture''. The technical side 
of the question concerns our ability to completely specify the geometry of a domain from the eigenvalues of 
its Laplacian. The question has been reinterpreted in the study of Schrodinger operators on metric graphs 
by Gutkin and SmUansky^ and restated in Algebraic Graphs Theory as "Which graphs are determined by their 
spectrum?" by van Dam and Haemers'*. (Through this work, the spectrum of a matrix M, denoted by Sm, is the set 
of its eigenvalues.) While we commonly employ different types of matrices to encode the structure of graphs, none 
has yet been shown to efficiently provide a complete graph invariant, i.e., a parameter that does not change under a 
permutation of the vertex labels. The spectrum of the adjacency matrix, for example, is a common invariant and 
easily seen to satisfy the "if" part of this statement; however, it is not a complete invariant, given the fact that co- 
spectral non-isomorphic graphs are abundant^' (see for instance Supplementary Information Section A). In the 
same spirit, physical scenarios have suggested various notions of refined spectra as a tool for distinguishing 
graphs, with partial degrees of success^"'. A common intersection for these approaches is Quantum Mechanics, 
arguably due to the popularization of quantum dynamics on graphs at the beginning of the last decade'". 

It is interesting, not only from the historical point of view, to observe that the strong link between Physics and 
graphs is via the Ising model, perhaps the most studied model in Statistical Mechanics. Originally proposed in 
1925" as a simplified description of the magnetic properties of materials, the Ising model has found a vast number 
of applications from Biology to Solid State Physics. Its great importance is emphasized by exact solutions and 
numerical techniques for the identification of phase transitions and critical phenomena''^. The Ising model 
framework seems particularly suitable to observe differences between Classical and Quantum Mechanics in terms 
of spectral information, since the quantum case is directly obtained by adding an appropriate (transverse) mag- 
netic field to the classical Hamiltonian. 

In what follows, we map a graph into an Ising model and interpret its energy spectrum as a graph invariant, 
before and after the "switch" from Classical to Quantum Mechanics. We demonstrate with exhaustive numerical 
examples that the quantum spectrum is a stronger invariant and propose a general framework to define physically 
meaningful graph polynomials. Determining whether the quantum energy spectrum is a complete invariant 
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remains an open problem. We perform experiments on a program- 
mable annealer with superconducting flux technology'\ Our pur- 
pose is to "hear the shape of an Ising model", by generating statistics 
of low energy states as the outcome of a noisy evolution. The experi- 
ment is run disregarding whether or not the state of the device 
follows an adiabatic path along its instantaneous ground state, there- 
fore against the prescription for successful annealing". In other 
words, we are not only interested in the ground state but, unconven- 
tionally, in the full output of a noisy computation. We obtain data on 
non-isomorphic graphs that are distinguished by their quantum 
energy spectra but not by the classical ones. 

Results 

Classical Cospectrality and Ising models. The Hamiltonian of the 
Ising modeP^ (or, equivalent, 2-state Potts model) on a graph G, with 
n vertices V(G) and edges £(G), is defined by the diagonal matrix 

H{GJ): = } H{iJ)^jY[MG)ljHii,j) (1) 

MeE(G) {i,)} 

where, for each edge {/, j}, H(i, j) ; = jH(fc) is a 2" X 2" matrix, 
with H(k) = (Tj, if fc = i, j and H{k) = I, otherwise. A{G) is the 
adjacency matrix with [A(G)],j = 1 if {i,j} e -E(G) and [A(G)]y = 
0, otherwise, is the Pauli matrix in the z-th coordinate axis, / is the 
identity matrix, and / is the strength of interaction. From now on, 
whenever the interaction strength is not expressly indicated as, e.g., 
in H{G), we implicitly set / = 1 for all edges. The partition function of 
the Ising model on G is 

Z(G,v)=Tr(e-/'«(°)), (2) 

where P'.= (fc^T}"' and v = e'*' — 1; kg is Boltzmann's constant, 
r e R"" is the temperature. By the Fortuin-Kasteleyn"" combina- 
torial identity, 2(G, v) is an evaluation of the Tutte polynomial"-^'', 
which is a fundamental invariant that determines many parameters 
including girth, chromatic number, etc. Remarkably, the Jones 
polynomial of a knot is contained in the Tutte polynomial'*'. Recall 
that, formally, G and G' are isomorphic if they are the same graph up 
to a relabeling of the vertices. This is denoted by G = G'. It is not hard 
to find graphs with the same Tutte polynomial (T-equivalent) that 
are not isomorphic^" "^': for example, all trees on the same number of 
vertices. 

Observe that two graphs G and G' have the same partition func- 
tion if and only if they share the same spectrum of the Hamiltonian in 
Eq. (1) (i.e. Z{G, v) =Z{G', v)<s>Sh(g) =Sh(g))- We say that G and 
G' are co-Ising if Sh(g) = Sh(g') • Since the Tutte polynomial is a gen- 
eralization of the partition function, if two graphs are T-equivalent 
then they share the same energy spectrum and thus are co-Ising. 
Thus, we know the following: 

G = G'=>Sh(g)=Sh{g'} 

(3) 

Sh(g) = Sh(g') f>G = G . 

Intuitively, we may attempt a refinement by adding a longitudinal 
field. The Hamiltonian of the Ising model on G with longitudinal field 
is defined by the diagonal matrix 

HL{G,J,h):=H{G,})+hM, (4) 

where M:=^"_^K(!) is a 2" X 2" matrix, with K{i) = 
®l^^Hi{k), Hdk) = a, a k = i and Hdk) = I otherwise. 
Physically hM can be interpreted as a constant external magnetic 
field applied to all vertices. Again, we set / = 1 and h = I unless they 
are explicitly indicated. We say that two graphs G and G' are longit- 
udinal field co-Ising if SHi^(Gj,h) —ShiIG'jJi) for values of / and h. 
The following equation summarizes what we know about graphs 



with this property (see Supplementary Information Section A and 
B for examples): 

C ^Zl'" c c ^ c 

^HiiGJ.h) - i>Hi{G'J,h)=^^H(GJ) - i>H{G'J) 

Shi(G) = Sh^(G')'^^Sh{G) = Sh(G') (5) 

W,^ SHi{G.J,h)=SHi{G'J,h)'^G=G'- 

From the diagonal matrices H{G) and M, we can define the energy 
and magnetization vectors as e„{G) = H{G)a_a and m^r = M„^at where 
(T = 0, 2, . . ., 2" — 1 runs over the classical states of the Ising model on 
G, where 0 denotes the ground state. With the use of these vectors, the 
bivariate Ising polynomial is defined as": 

Z(G;x,;.)=^x^"(°V-". (6) 

Notice that the spectrum Sfji^(Gj,h) can be obtained from Z(G; x,y) for 
all values of the constants / and h, since a change in these parameters 
is just a rescaling of the coefficients x and y. The Ising polynomial 
generalizes the partition function in Eq. (2) because Z{G, e"''', 1) = 
Z{G, e"' — 1), encodes the matching polynomial, is related to the van 
der Waerden polynomial, and is contained in a more general poly- 
nomial introduced by Goldberg, Jerrum and Paterson^''"^^. The 
bivariate Ising polynomial in Eq. (6) can be intuitively generalized 
by working with any physical observable in addition to energy and 
magnetization. If we denote by the eigenvalues of a diagonal 
matrix (or observable) A', we can then define a multivariate poly- 
nomial 

Z{G; X, y, z,) = ^ x^'{G)y-' Yl zf . (7) 

a 

An example is given by the (permutationally invariant) spin-glass 
order parameter used by Hen and Young"^*". 

Quantum Cospectrality. The invariants that we have so far 
considered belong to Classical Physics. We can now move into a 
quantum mechanical regime by adding a further field. The 
Hamiltonian of the quantum Ising model on G, as proposed by 
Lieb, Schultz, and Mattis^' (see also'"') is defined by the matrix 

HT{G,Lh, A) ■. = HL{G,I,h)-\-AMr, (8) 

where A e IR is a transverse external magnetic field; here Mj : = 
J21^i T{i) is a 2" X 2" matrbc, with T[i) = ^l^-^Hjik), Hjik) = 
a^ifk = i and Hjik) = I otherwise. As in the longitudinal case, Mj 
does not depend on G. Two graphs G and G' are said to be quantum 
co-Ising if SHy(G,;,fi,A) =Sht(g',/,a,A) for all values of /, h and A. It 
follows from the definition that two graphs are quantum co-Ising 
if they are isomorphic. The quantum partition function is defined 
analogously to the classical one: 

Zt{G, /i, 7, h, A)=rr(e-I'"^^'^'''''''^^y (9) 

Two graphs are quantum co-Ising if and only if they have the same 
quantum partition function. The "if" part of this statement comes 
directly from the definition. For the "only if part, observe that in the 
limit p cc, Zt{G) ^ X'oe^^^" determines the lowest eigenvalue Eq 
with its multiplicity Vq. Similarly, in the same limit Zr(G)e'*^°/ 
determines the value and multiplicity of the second smallest 
eigenvalue. The whole spectrum is obtained iteratively. The 
statement above and its proof are valid only for systems of finite 
size. It is a well-known fact that different Hamiltonians can have 
the same partition function in the thermodynamic limit. 

We tested numerically the converse of this fact by computing 
the smallest eigenvalue for h = J = A = 1. We tested all graphs with 
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n £ 9, all bipartite graphs with « £ 11, all vertex transitive graphs 
with « £ 15, all regular graphs with n £ 1 1, and all trees with n £ 14 
(also considered in^^). We failed to find a counterexample. Hence, 

(10) 

Sht{g) =Sht{g') =>G = G'. 

The transverse field Ising HamUtonian is a sum of non-commuting 
terms and determining its full spectrum requires the diagonalization 
of a 2" X 2" matrix. We cannot generalize directly the quantum 
partition function to a generating Ising polynomial as done in the 
classical case - when eigenvalues are integers (for J = h = 1) - 
although we can use the well-known Suzuki-Trotter formalism to 
obtain a classical approximation^'; the direct calculation of the eigen- 
values is notoriously expensive, due to the size of the problem, and 
prone to errors, making it difficult to numerically show the existence 
of non-isomorphic quantum co-Ising graphs. A reasonable first 
approximation for this task is to compute the absolute largest eigen- 
value. That is what we have done in our tests. Taking into account 
such difficulties, finding non-isomorphic quantum co-Ising graphs 
is an open problem. Natural candidates are graphs for which iso- 
morphism testing is known to be harder to solve {e.g., graphs 
for which the Weisfeiler-Lehman algorithm fails)'". Nevertheless, 
we emphasize that spectral information provided by Quantum 
Mechanics is more accurate than Classical Mechanics. It must be 
said that there are only a few precise (and in fact negative) statements 
about the physically inspired graph invariants which have been 
introduced recently (see' and the references therein) and that 
purely numerical analysis does not guarantee sufficient generality. 

Experiments. Disregarding computational complexity aspects, we 
have highlighted that from the theoretical point of view one can 
hear the shape of certain quantum Ising models, while it is not pos- 
sible for the classical analogue. We subsequently encode on the same 
physical system pairs of non-isomorphic graphs that are co-spectral, 
longitudinal field co-Ising (and consequently co-Ising), but not 
quantum co-Ising. Rather remarkably, our set up finds an experi- 
mental implementation in the optimization technique called quantum 
annealing^*-^^'^^ . In this technique, the system evolves adiabaticaUy 
according to the following time-dependent Hamiltonian 

Hqa{GJ, h. A, s)=sHi{G,J, ft) + (l-s)AMr, (11) 

where s = t/Ttot, t is a time parameter and Tfot is the total duration of 
the dynamics. At the beginning of the computation, the system is 
prepared in the ground state of the initial simple Hamiltonian 
Hqa^G, J, h. A, 0) = AMj-. On the basis of the adiabatic 
theorem'*, adiabatic quantum annealing with general HamUtonians 
has been shown to be a universal model of computation by Aharonov 
et aV^. In synthesis, the core idea is to evolve the system slowly 
enough towards a final ground state, which is the solution of a 
computational task. While the success of this paradigm depends 
on the ability of avoiding level crossings with ad hoc annealing 
schedules, Brooke et aU^ experimentally observed that tunneling 
can hasten convergence to the solution. 

In the setting specified by Eq. (11), we are interested in measuring 
the observables eo, mg, and Oq. In a realistic situation, temperature 
and environment will usually excite the system. While these effects 
are disruptive in the standard applications of quantum annealing, we 
regard such a non-ideal implementation as a way to generate the 
statistics of low energy states on which we measure the correspond- 
ing observables. For this purpose, we run experiments on a D-Wave 
Vesuvius programmable annealer. The hardware consists of 503 
usable logical bits on an integrated circuit with superconducting flux 
qubits (see" for details on the technology). Quantum effects on the 
chip are currently under investigation and there is evidence of 
quantum annealing on random spin glass problems'" The 



Hamiltonians that can be realized with the device are exactly of the 
type in Eq. (11), where s is a non-linear function of time. The most 
general form of the final Hamiltonian Hi is given by an Ising model 
whose possible spin interactions are constrained by the chip archi- 
tecture. A particular limitation of the hardware is that measurements 
can be performed only at the end of the evolution. Thus, the maximal 
information that we can extract is encoded in the multivariate poly- 
nomials of Eq. (7). On the other hand, the final state of the chip is a 
result of a dynamics also governed by the transversal field Mj. In fact, 
our experiments attempt to identify the effects of Mj- in the final 
statistics after the measurement outcomes are filtered out using vari- 
ous type of multivariate polynomials. 

We have tested the annealer on two pairs of non-isomorphic 
graphs G and G' (G13 and G'n, G27 and G'27, in Supplementary 
Information Section B and C) such that Sa(g) = Sa(g')' Shl^gjm) ~ 
Shi(g',/,a). and Shj{g) ^Sh,(g')' '■^■> with equal spectra of the adja- 
cency matrix, equal classical spectra, even with a longitudinal field, 
and different quantum spectra. To illustrate a possible (arbitrary) 
refinement as introduced by Eq. (7), we include an extra observable, 
SV, corresponding to the next-nearest neighbor interaction energy: 
n.''=^,,H{i,j)[A''{G)]... Notice that H(G) = nV2. Figure 1 

shows the statistics of measurement outcomes when the states are 
distinguished through the doublet of observables {energy, magnet- 
ization} on the pairs {G13, G'13}, for J = h = 1/7 and J = h = 1. These 
are respectively the smallest and the largest values that can be reliably 
set on the hardware. The final states are organized according to the 
values of the two observables. As a consequence of the fact that the 
two graphs are co-Ising, the measured values of the pairs {energy, 
magnetization} are the same, and cannot be used to distinguish the 
two graphs. Moreover, the shape of the two distributions is also the 
same up to statistical errors. The shape of this distribution is assumed 
to be a consequence of (noisy) open system quantum dynamics''"'" 
(see Supplementary Information Section D for a comparison bet- 
ween experimental and thermal statistics). This means that we are 
not able to identify differences in the final distributions that may 
arise due to the different quantum spectra, i.e. due to non-equivalent 
quantum evolution along the annealing schedule. 

The graphs are indistinguishable by measuring energy and mag- 
netization only. However, they become distinguishable in Figure 2 by 
measuring the triplet {energy, magnetization, SV}, as clearly visible in 
the statistics obtained with the chip. The pair {G27, G'27} is not dis- 
tinguished by the triplet on the experimental data, as showed in 
Figure 3. It should be possible, in principle, to classically distinguish 
these graphs with the introduction of additional observables. 
Similarly to what happens for the G13 pair, there are no noticeable 
differences in the shape of the final distributions that can detect 
differences in the quantum spectra. 

Conclusions 

The interplay between combinatorics and the classical Ising model is 
well-established. We have introduced a general family of physically 
meaningful graph polynomials suggesting a hierarchy of graph 
invariants. We have demonstrated that the quantum Ising model is 
a finer sieve to distinguish graphs than its classical analogue by con- 
sidering the quantum partition function as a graph invariant. We 
have tested experimentally its distinguishabUity power on a D-Wave 
programmable annealer, by taking graphs with different quantum 
spectra and the same classical Ising partition function. We used the 
hardware unconventionally to generate the statistics of low energy 
eigenvalues rather than focusing on the ground state. The data 
obtained can distinguish one pair of graphs when measuring with 
respect to a classical refinement of the partition function. We did not 
find any measurable difference in the statistics of measurement out- 
comes of the two pairs that can be related to non- equivalent quantum 
dynamics. Notice that the transverse field spectra are very similar 
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Figure 1 | The statistical distribution of measurement outcomes JV^^^ on the pairs {613,613} obtained by averaging over 100 cycles for each of the 100 
different embeddings considered (10000 programming cycles in total). The horizontal red line corresponds to the median of the data while the edges of 
the blue boxes correspond to the 1st and 3rd quartile. Each red cross is an outlier measurement. The outcomes have been filtered after choosing the pair of 
classical observables {eo, mo}. Data showed in the left panels correspond to the choice J = h = 117. In the right panels / = h = 1, that is the maximum 
strength of the couplings allowed by the hardware. With the given choice of classical observables, the distribution of measurement outcomes is not able to 
distinguish the two graphs, nor at the classical, neither at the quantum level. 



(Fig. 9 in the Supplementary Material). Of course, differences ex- 
pected in an ideal quantum system are possibly lost due to decoher- 
ence when approaching the classical regime at the end of the adiabatic 
evolution. 

Going beyond the scope of this work, it would be interesting to 
compare the experimental data with numerical simulations of the 
corresponding open quantum spin system at finite temperature*'. 
We propose two approaches to amplify the differences in the quan- 
tum spectra: (a) reduce substantially the annealing time; (b) perform 
measurements when the transverse field is on. Both approaches re- 
quire a modification of the current control of the hardware. Another 



interesting goal is to define other efficient observables, such as Q^, 
that would amplify the possible differences in the measurement stat- 
istics. From the theoretical point of view a natural open question is 
whether the transverse field alone is sufficient to define a complete 
spectral graph invariant. 

Methods 

Experimental data collection. In order to collect enough statistics for averaging over 
biases and systematic errors, we have considered 100 embeddings in the chip for each 
graph. To average over precision errors when setting the intended couplings on the 
machine, we have run 100 programming cycles for each embedding. For each cycle, 
we have performed 1000 measurements. All the experiments have been performed 
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Figure 2 | The statistical distribution of measurement outcomes N^cc on the pairs {613,613}. The outcomes have been now filtered after choosing a 
triplet of classical observables {eo, mo. Or]. Data showed in the left panels correspond to the choice J= h = 111. In the right panels ] = h = 1 . Using a third 
observable distinguishes the two graphs at the classical level. 
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Figure 3 | The statistical distribution of measurement outcomes A/q^.^ for the pairs {G27, G'27}./— h— 1/7 in the left panels. /— h— 1 in the right panels. 
Filtering the outcomes after choosing the triplet {cq, mo, O.^} does not distinguishes the graphs at the classical level, and the introduction of additional 
observables is needed. The shape of the two distributions is also the same, meaning that the two graphs are not distinguished at the quantum level either. 



choosing the shortest annealing time allowed by the hardware {Ttat — 20 fis) in order 
to minimize the effects of thermal excitations. 
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